Testing the self-duality of topological lumps in SU(3) lattice gauge theory 
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We discuss a simple formula which connects the field-strength tensor to a spectral sum over certain 
quadratic forms of the eigenvectors of the lattice Dirac operator. We analyze these terms for the 
near zero-modes and find that they give rise to contributions which are essentially either self-dual 
or anti self-dual. Modes with larger eigenvalues in the bulk of the spectrum are more dominated by 
quantum fiuctuations and are less (anti) self-dual. In the high temperature phase of QCD we find 
considerably reduced (anti) self-duality for the modes near the edge of the spectral gap. 
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Chiral symmetry breaking is one of the central pil- 
lars in our understanding of QCD. For the underlying 
mechanism an interesting picture based on instantons 
has been developed during the last 20 years [Q. The 
relevant excitations in the QCD vacuum are believed to 
consist of a weakly interacting fluid of instantons and 
anti-instantons. For a single instanton or anti-instanton 
the Dirac operator has an exact zero-mode. For the 
weakly interacting instantons and anti-instantons instead 
of many zero modes one expects the build-up of a non- 
vanishing density of small but non-zero eigenvalues cor- 
responding to so-called near zero-modes. This spectral 
density near the origin is related to a non-vanishing chi- 
ral condensate through the Banks-Casher formula |^]. 

During the last year we have seen several papers 
1^ which test the scenario of chiral symmetry breaking 
through instantons with ab-initio calculations on the lat- 
tice, some supporting the instanton picture, some chal- 
lenging it. The common basic idea is to study the prop- 
erties of the near zero-modes of the lattice Dirac oper- 
ator and to analyze whether one can identify localized 
structures which resemble instantons. The underlying 
assumption is that the near zero-modes are only slight 
deformations of zero-modes which in turn are known to 
be localized at the same position as the instanton. 

In this letter we derive a simple formula which relates 
the field strength tensor Ff^^, to a spectral sum over cer- 
tain quadratic forms of the eigenvectors of the lattice 
Dirac operator. The terms in this spectral representa- 
tion of the field strength, i.e. the quadratic forms of the 
eigenvectors, have a clear interpretation. One can ana- 
lyze the contributions of the near zero-modes and test the 
scenario of chiral symmetry breaking through instantons. 

In particular we will focus on analyzing the duality 
properties of the field strength. The possibility that the 
localized lumps in the field strength are not (anti) self- 
dual, i.e. do not have a fundamental property of (anti) 
instantons, has been brought up in Ref. , a paper which 
challenges the local quantization of topological charge in 
integer units. However, this possibility is ruled out by 
our current results: We find that the contributions to Ff^„ 
which come from the near zero-modes do indeed show a 



high degree of (anti) self-duality. For contributions from 
eigenmodes with larger eigenvalues we find a decrease of 
(anti) self-duality as they are more dominated by quan- 
tum fiuctuations. 

When the temperature is increased beyond its criti- 
cal value, QCD undergoes a phase transition where chi- 
ral symmetry is restored. In the instanton picture this 
transition is believed to be related to the formation of 
tightly bound instanton anti-instanton molecules. The 
spectrum develops a gap and due to the Banks-Casher 
formula the chiral condensate vanishes. When analyz- 
ing the modes with eigenvalues close to the edge of the 
spectral gap, we do indeed find a considerably reduced 
amount of (anti) self-duality as is expected for the tightly 
bound molecules. 

We will denote the Dirac operator hy D — 7/^-0^ with 
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p. The gauge potential is anti-hermitian 
and in components is given by A^^ = —iA^jTa where 
Ta, a = 1..8 are the generators of su(3). In this notation 
the field strength tensor is given by F^^ = [D^, D^], and 
in components reads = —iF^^^Ta- When evaluating 
the square of the Dirac operator one finds 



(1) 



The field strength F^^, can be projected out by multiply- 
ing with 

7/i7iy Bind taking the trace over the Dirac 
indices. This results in the formula (/i 
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for the field strength F^^. For a lattice Dirac operator 
the situation is slightly more involved. On the lattice the 
Dirac operator is a difference operator and has two space- 
time indices x, y. When repeating the above calculation 
for e.g. Wilson's lattice Dirac operator Dw, one finds 
(the indices fi,iy on the right hand side are not summed) 
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Thus one finds that the field strength is slightly smeared 
out among several lattice points. We remark that an 
equivalent formula holds for staggered fermions where 
liilv is replaced by a staggered factor depending on /i, v. 
For a general Dirac operator, such as the overlap operator 
0, the fixed point operator |^ or the chirally improved 
operator |^ (which we use here) one finds (from now on 
we drop "^(a^)") 
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with some function ty{x)f^i, which describes the smear- 
ing out of the field strength around the central point x. 
For the overlap operator this function is fast decaying 
in I a; — 2/1 but has infinite support while for a finite ap- 
proximation of the fixed point operator or the chirally 
improved operator ^^(a;)^^ will vanish for a sufficiently 
large separation of x and y. One can show that a sub- 
stantial contribution comes from x — y and below we will 
concentrate on this particular case. 

We now use the spectral representation of the Dirac op- 
erator to express the right hand side of Eq. (||) in terms 
of the eigenvectors D{x,y) = J2j y{j\- In this 

formula we denote the eigenvectors by \ j)x and write ex- 
plicitly only the space-time index x, while the color and 
Dirac indices are denoted with the bra-ket notation. The 
corresponding eigenvalues are denoted by . We remark 
that for non- normal operators, the bra's y(j\ have to be 
replaced by left eigenvectors and only the \j)x are the 
conventional right eigenvectors 

After inserting the spectral decomposition of D into 
our formula (^) we find: 



-2 !/0l7M7.r,|j% . (6) 

In order to project onto the a-th component of the field 
strength tensor we have used = i2 TrcfF^i/Ta]. 

Equations (||) and (|^) provide the announced spectral 
decomposition of the field strength tensor. One can now 
study the individual contributions f^iyix)y,j in order to 
analyze properties of the field strength tensor. It is in- 
teresting to note that the exact zero modes do not con- 
tribute to the spectral decomposition (||),(^. Our for- 
mula has been tested on lattice discretizations of contin- 
uum instantons 

Let us briefly comment on the setting of our numerical 
calculation. We study QCD in the quenched approxima- 
tion using the Liischer-Weisz action with coefficients 
from tadpole improved perturbation theory. We work 
on lattices of size 16^ for the zero temperature ensem- 
bles and on 6 X 20^ lattices for the ensembles in the high 
temperature chirally symmetric phase. The leading cou- 
pling /3i of the Liischer-Weisz action is (3i = 8.45 which 
corresponds to a lattice spacing of a = 0.095 fm 




FIG. 1. Contributions fi2 (top plot) and f^^ (bottom 
plot) of the first near zero-mode to the field strength. We 
show a slice of a 16* lattice. 



This value gives a temperature of 346 MeV for the high 
temperature ensemble. We use the chirally improved op- 
erator which is a systematic expansion of a solution of 
the Ginsparg- Wilson equation The computation of 
the eigenvalues and eigenvectors of the Dirac operator 
was done with the implicitly restarted Arnoldi method 

Let us begin our analysis of the terms in the spec- 
tral representation (^,(||) by displaying a single contri- 
bution f^^{x)yj and compare it to its dual f^,j{x)yj = 
1/2 £^ivpafpa{^)y,j- In particular in Fig. ^ we show a 
slice through the 16^ lattice and display the /fj compo- 
nent in the top plot and the component in the bottom 
plot. We have set x — y, i.e. we show the center of the 
smeared-out field strength, and we evaluate the contri- 
bution for the first near zero-mode, i.e. the eigenvector 
with the smallest non-vanishing eigenvalue. 

Both plots show a large peak accompanied by smaller 
wiggles. The peak changes sign but otherwise essentially 
keeps its size and shape when going from the to the 
/34 component. Thus the peak is anti self-dual. For 
the smaller wiggles it is not straightforward to identify 
such a simple duality behavior - they are mainly quan- 
tum fluctuations. When inspecting many of these plots 
we find that the typical pattern consists of small quan- 
tum fluctuations and pronounced peaks which are either 
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^IG. 2. Tr Ff,^{xyF^^{xY (top plot) and Tr F^^{x)' 
Fij,v{x)' (bottom plot). The primes indicate the that we 
take into account only the contributions of the lowest 6 near 
zero-modes as defined in Eq.(^. 



self-dual or anti self-dual. For contributions of eigenvec- 
tors with larger eigenvalues, so called bulk-modes, one 
finds that the abundance of isolated peaks decreases and 
the contributions become more dominated by quantum 
fluctuations. 

Quantities of more physical interest are the action den- 
sity Tr FpjyFp^ and the topological charge density Tr 
Ff_ii,F^,^. In particular a direct comparison of the two 
quantities shows if the field is (anti) self-dual in all com- 
ponents. Again it is most interesting to analyze the con- 
tributions of the near zero-modes and see whether they 
are dominated by (anti) self-dual lumps as expected from 
the instanton picture. To this purpose we truncate the 
sum over the eigenvalues Xi in Eq. (|^) and take into ac- 
count only the first TV eigenvalues. In Fig. || we show 



TrF^,(x)'F,,,(x)' 
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in the top plot. The color indices a = 1..8 on the right 
hand side of (Q) are summed to produce the color trace 
on the left hand side. The corresponding contribution to 
Tr F^v{x)' F^ij{x)' is constructed in the same way and is 
displayed in the bottom plot of Fig. ^. For the figure we 
set TV = 6, i.e. we take into account only the contributions 
from the lowest 6 near zero-modes. 

Both plots in Fig. ^ show two pronounced peaks 
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with one of them changing sign when going from Tr 
'^uixYF^.^ix)' to Tr F^^{x)'F^^{x)' . We remark that 
igs. |l| and|| were made from the same gauge configura- 
tion and the anti self-dual sign-changing peak of Fig. ^ 
is the same peak which we already saw in the contribu- 
tion to the 5-component displayed in Fig. |l|. The other 
smaller peak in Fig. |2| is a self-dual fluctuation in a dif- 
ferent color component. Since the composed quantities 
Tr F^^ix)' Ff,^{x)' and Tr F^^{xyF^^{xy are built from 
products of the single components the relative size 
of the quantum fluctuations is suppressed considerably 
compared to the larger, (anti) self-dual structures. We 
find that the contributions of thejiear zero-modes to Tr 
F^_i,v[xy F^^y{xy and Tr F^,y(a:)'F^y (x)' are entirely dom- 
inated by lumps which are either self-dual or anti self- 
dual. 

In order to go beyond an illustration of the duality 
properties by a few examples we now discuss an observ- 
able which allows to test (anti) self-duality systemati- 
cally. Its construction is similar to the chirality observ- 
able proposed in M. We define the ratio 



r{x) 



Ti F^,{xy F^,{xy - Ti F^,{xy F^,[xy 



(8) 



TTF^,{xyF^,{xy + TvF^,{xyF^,{xy 

For a space-time point x where the gauge field is self- 
dual the numerator will vanish while the denominator is 
finite and r{x) equals to 0. Conversely for an x where the 
gauge field is anti self-dual the role of numerator and de- 
nominator are exchanged and r{x) = oo. For space-time 
points without definite duality properties r{x) assumes 
some finite value between and oo. The transformation 
R{x) — 4/7rarctanr(a;) — 1 maps the interval [0, oo) into 
the interval [—1, 1]. For configurations which are domi- 
nated by (anti) self-dual lumps one expects values near 
±1. As for the local chirality variable of jj] it is inter- 
esting to study different selections for the lattice points 
X in R{x). In Fig. ^ we use all lattice points (top curve), 
the subset of 50% of the lattice points supporting the 
highest peaks of |Tr F^^{xy F^j^v{xy\ (middle curve) and 
also a cut of 10% (bottom curve). Again we use the 6 
lowest modes in the series for F^^(xy . The histograms 
were computed by averaging over 50 configurations. 

Fig. ^ shows that the contributions of the near zero- 
modes to the spectral representation of the field strength 
are highly (anti) self-dual. Even when including all lat- 
tice points (100%, top curve) one finds a pronounced dou- 
ble peak. When throwing away 50% of the lattice points 
with small |TrF^,y(a;)'F^y (a;)'| this cuts mainly into the 
center of the histogram, showing that the most (anti) self- 
dual excitations are indeed the large peaks sticking out 
of the quantum fluctuations. This trend continues when 
focusing on only the highest 10% of the peaks. Since our 
analysis is based on the exact spectral decomposition of 
the fleld strength it truly reflects self-duality of the topo- 
logical lumps in the gauge field. 

It is interesting to perform the same analysis also in 
the chirally symmetric phase of QCD. As discussed, the 
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FIG. 3. Histograms for the observable R which tests the 
(anti) self-duality of the near zero-modes (compare Eq. (^)). 
We display our data without a cut (top curve), a cut of 50% 
(middle curve) and a cut of 10 % on the number of lattice 
points supporting the highest peaks of \TrFi_i^{x)' F^^{x)'\. 
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FIG. 4. The same as in Fig. ^ but now for the ensemble in 
the chirally symmetric phase of QCD. 



instanton picture describes this phase by tightly bound 
molecules of instantons and anti-instantons |l[ . In Fig. ^ 
we show our results for the histograms of the duality 
observable R for our ensembles in the high temperature 
phase. Again we average over 50 configurations. F^^{xy 
is now computed from the 6 modes with eigenvalues clos- 
est to the edge of the spectral gap. Obviously the amount 
of (anti) self-duality has decreased considerably which is 
exactly what is expected to happen when instantons and 
anti-instantons condense to tightly bound molecules and 
lose their identity as isolated (anti) self-dual objects. 

Let us briefly summarize the obtained results. We de- 
rive a representation of the field strength in terms of a 
spectral sum of quadratic forms of the eigenvectors of the 
lattice Dirac operator. The contributions of individual 
eigenvectors to the field strength can be analyzed. We 
find that the near zero-modes are dominated by (anti) 
self-dual lumps, while modes with larger eigenvalues fur- 
ther in the bulk of the spectrum become dominated by 
quantum fluctuations. Our spectral representation thus 
decomposes the field strength tensor into contributions 
(the near zero-modes) sensitive to the large scale fluctu- 
ations and contributions (the bulk modes) dominated by 
quantum fluctuations. In the high temperature phase we 
find a reduced amount of (anti) self-duality for the modes 



with eigenvalues near the edge of the spectral gap. Our 
results support the picture of chiral symmetry breaking 
by instantons and its restoration through the formation 
of instanton anti-instanton molecules. 
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